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1.  Introduction: 


Subject  report  was  prepared  by  Measurement  Analysis  Corporation 
(MAC)  under  contract  for  6585th  Test  Group.  Its  principal  purpose  was 
the  design  of  a  high  precision  test  platform  suitable  for  the  testing 
of  future  generations  of  inertial  guidance  components  and  systems.  This 
platform  also  will  serve  as  a  prototype  facility  for  test  concept 
verification  for  a  planned  Precision  Guidance  Test  Facility. 

This  evaluation  of  the  report  was  performed  by  Or.  Martin  G.  Jaenke 
serving  as  a  consultant  in  a  temporary  employment  limited  to  200  hours. 
The  emphasis  of  the  evaluation  was  on  predicting  the  performance  of  the 
proposed  system  under  realistic  assumptions  of  seismic  disturbance  inputs 
and  on  the  critical  analysis  of  some  design  details  which  are  not 
covered  sufficiently  in  subject  report. 

Due  to  the  limited  time  available,  all  numerical  calculations 
required  for  the  analysis  were  performed  by  using  manual  methods 
(pocket  calculator)  because  it  was  felt  that  the  development  of  new  or 
the  adaptation  of  existing  computer  programs  would  be  too  time  cr  su- 
ming.  This  involved  certain  simplifications,  however,  the  results 
obtained  provide  at  least  an  "order  of  magnitude"  picture  of  trends. 

In  the  following  paragraphs  the  notation  p.  XX  refers  to  page  XX  of 
subject  report.  To  facilitate  typing,  all  sketches,  formulas  and 
tables  are  comprised  on  "plates",  refered  to  a  FI  YY. 

Appendix  A  contains  a  short  discussion  of  the  possible  approaches 
to  present  variables  which  are  only  statistically  defined,  a  problem 
which  frequently  led  to  confusion. 

Appendix  B  is  a  copy  of  a  paper  which  the  author  presented  to  the 
AIAA  in  197^.  Besides  being  useful  for  explaining  the  assumptions  of 
seismic  disturbances  made  in  the  evaluation,  it  is  thought  to  be  valu¬ 
able  in  the  planning  process  for  a  future  Precision  Guidance  Test 


The  system  is  described  in  subject  report  on  p.  3  through  p.  46. 

PI  1  gives  a  short  description  sufficient  for  the  understanding  of  the 
discussion  in  the  following  paragraphs.  It  shows  that  the  system  is  a 
dual  isolation  system.  The  primary  reaction  mass  M2  carries  the  pay- 
load  PL  (test  stand  and  test  item)  and  the  instrumentation  I  necessary 
to  sense  the  motion  of  the  mass.  The  secondary  reaction  mass  Ml  is  a 
necessary  prerequisite  to  make  the  dual  isolation  concept  possible. 

The  advantage  of  this  concept  is  that  it  provides  significantly  higher 
high  frequency  attenuation  of  seismic  disturbance  inputs  from  the  ground 
than  a  simple  mass-spring  system.  There  are  two  separate  sets  of 
springs,  S.  The  spring  action  is  provided  by  pneumatic  isolators  (Barry 
Serva-Level  Pneumatic  Isolators  #133-12).  The  vertical  (S y)  and  hori¬ 
zontal  (S^)  spring  actions  are  generated  within  the  same  unit.  The  ver¬ 
tical  action  is  truly  pneumatic  while  the  horizontal  action  is  due  to 
the  elastic  properties  of  a  diaphragm  which  supports  the  load  carrying 
element  (for  detailed  description  see  p.  C-8  and  p.  K-34) .  Electomag- 
netic  actuators  A  operate  on  the  primary  reaction  mass  following 
commands  derived  from  the  instrumentation  I  and  processed  by  appro¬ 
priate  servo  loop  elements. 

The  performance  of  the  proposed  platform,  to  the  extent  that  it  is 
relevant  for  the  foliowing  discussion,  is  described  in  PI  2  -  PI  4. 

This  performance  description  is  in  terms  of  the  transmi ss i b i 1 i ty  between 
the  disturbanoe  motion  input  from  the  ground  and  the  resulting  motions 
of  the  primary  reaction  mass.  "ine  information  is  taken  directly  from 
subject  report  and  is  accepted  at  face-value  at  this  point  of  discus¬ 
sion.  A  proposed  revision  of  the  performance  indicators  is  described 
in  para  3.3  of  this  evaluation  and  its  consequences  for  testing 
performance  will  be  discussed. 
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PI  2  and  3  describe  the  t  ransmi  ss  ib  i  1  i  ty  for  translational 
motions  in  the  horizontal  and  vertical  directions.  PI  2  shows  first  as 
a  function  of  frequency  the  "passive"  transmissibi 1 i ty ,  i.e.  the  per¬ 
formance  of  the  system  with  the  active  control  servo  loop  not  in 
operation.  Next,  both  show  the  required  transmiss ibi 1 i ty  which  is 

necessary  to  reduce  the  disturbance  inputs  which  were  actually  observed 

-8 

in  AITL  to  the  specified  level  of  10  g.  These  disturbance  inputs 
agree  reasonably  well  with  observations  made  previously  by  the  author 
and  derived  from  the  literature.  Finally,  PI  2  and  3  present  the 
transmi ss ib i 1 i ty  predicted  with  the  active  stabilization  system  in 
operation.  It  is  shown  that  the  system  specif ications  with  respect  to 
translat ional  motions  are  expected  to  be  met,  a  result  which  is 
intuitively  accepted  by  this  author.  This  is  not  the  case  with  respect 
to  rotational  (angular)  motions  which  have  signi f icantly  more  severe 
consequences  on  gyro  and  accelerometer  tests  than  translational  distur¬ 
bances.  Unfortunately,  the  whole  area  of  rotational  transmi ss ib i 1 i ty 
is  treated  only  superfically  in  subject  report  and  therefore  will  be 
the  object  of  a  critical  analysis  in  para  3-3  of  this  evaluation. 

PI  k  presents  the  expected  angular  transmi ss ib i 1 i ty  from  subject 
report.  It  will  be  noted  that  the  plate  does  not  contain  information 
on  passive  and  required  t ransmi ss i b i 1 i ty .  The  former  can  not  be  found 
anywhere  in  subject  report  and  the  latter  deficiency  is  due  to  the  fact 
that  the  contractor  did  not  use  information  on  angular  disturbance  in¬ 
puts  although  it  was  available  to  him.  However,  as  stated  above,  PI  4 
will  be  used  at  face-value  for  the  following  analysis  of  the  effects 
of  angular  motions  on  test  accuracy.  PI  5  shows  the  angular  distur¬ 
bance  inputs  in  terms  of  a  power  density  spectrum  which  must  be  expected 
for  the  operation  of  the  proposed  platform.  It  is  taken  from  Fig  5  of 
Appendix  B  where  detailed  information  on  its  conceptional  background  can 
be  found.  The  essential  features  of  this  spectrum  are  repeated  as 
curve  A  in  PI  b,  but  here  in  terms  of  an  amplitude  spectrum  (for  a 
discussion  of  these  different  types  of  presentation  see  Appendix  A). 
Curve  B  of  PI  6  shows  the  resulting  platform  motions  obtained  after 
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attenuation  by  the  active  isolation  system  in  accordance  with  PI  4. 

P)  7  shows  the  growth  of  angular  platform  deflection  with  time  derived 
from  the  spectra  of  PI  6,  curve  A  without  and  curve  B  with  the  active 
isolation  system  in  operation.  It  is  of  interest  to  note  that  B  attains 
a  final  value  of  about  0.02  sec,  a  value  which  served  as  a  specifi¬ 
cation  for  the  design  contract.  PI  8  shows  the  spectra  of  drift  induced 
by  platform  angular  motions  into  a  gyro  under  test,  again  without  and 
with  platform  stabilization.  There  are  two  mechanisms  causing  this 
drift:  one,  essentially  inversely  proportional  to  frequency  at  low 
frequencies  is  the  "pick-up"  of  components  of  the  earth  rotation  vector 
due  to  angular  deflections  of  the  test  item.  The  other,  essentially 
constant  at  higher  frequencies,  is  due  to  the  sensing  of  angular  rates 
by  the  input  axis  of  the  gyro  under  test.  The  total  rms  drift  value  in 
the  frequency  band  from  10  to  30  Hz  is  4.3  x  10  3  deg/hr  without  and 
5  x  10  **  deg/hr  with  active  isolation.  These  values  are  obtained  by 
numerical  integration  of  the  power  spectra  which  correspond  to  the 
amplitude  spectra  of  PI  8.  It  is  seen  that  even  with  active  isolation 
the  platform  does  not  reach  a  goal  of  10  ^  deg/hr  a  goal  which  was 
obtained  from  PI  9  which  is  taken  from  Ref  1. 

A  similar  situation  prevails  in  accelerometer  testing.  PI  10 

shows  the  accelerations  induced  by  angular  platform  deflections  which 

lead  to  the  "pick-up"  of  components  of  the  gravity  vector.  The  total 

rms  value,  obtained  in  the  same  way  as  with  PI  8,  while  improving  from 

9  x  10  ^  to  3  x  10  ^  Sons  by  using  the  active  isolation,  does  not  reach 
-8 

a  goal  of  10  9rms  which  was  obtained  from  PI  11,  again  taken  from 

Ref  1.  It  is  important  to  note  that  the  acceleration  errors  induced 

-8 

by  angular  platform  motion  can  not  be  reduced  to  the  goal  of  10  g 
while  the  errors  induced  by  translational  platform  motions  apparently 
can,  as  was  discussed  in  connection  with  PI  2  and  PI  3  above. 

While  the  proposed  platform  does  not  reach  the  ideal  goals  formul¬ 
ated  in  connection  with  the  testing  of  very  advanced  future  inertial 


RMS  DRIFT  IN  BOND  FROM 
10-8  HZ  TO  30  HZ 
A  :  4,3  x  10  '3  Lpeg/hrJ 

B  :  5  X  10"5  G>EG/wJ 
GOAL:  10  ^  (deg/i^RM! 
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instruments,  it  is  of  interest  to  compare  its  performance  with  that  of 
precision  platforms  either  in  existence  or  under  development  at  other- 
test  facilities.  PI  12  provides  this  comparative  information,  curve  A 
for  the  proposed  platform  (taken  from  PI  A)  and  curve  B  for  the  existing 
platform  at  the  Frank  J.  Seiler  Research  Laboratory  (information  taken 
from  Fig.  6  of  Ref.  2  or  Fig.  7  of  Ref.  3).  Curve  C  presents  the  results 
obtained  from  a  prototype  precision  platform  which  was  developed  for 
the  now  discontinued  NASA  Electronic  Center  in  Cambridge,  MA.  The 
information  is  taken  from  Fig.  31  of  Ref.  A.  A  and  B  are  equal  in  the 
high  frequency  range  but  B  is  inferior  in  the  middle  range  which  con¬ 
tributes  quite  significantly  to  the  overall  rms  test  error  value.  C  is 
inferior  to  both  except  in  the  low  frequency  range.  Thus,  is  seems 
that  the  proposed  system  is  at  least  compatible  with  the  apparent  state 
of  the  art  and  superior  to  other  systems. 

3.  Detail  Analysis: 

3.1  Induced  Rocking  Motions 

The  structural  lay-out  of  the  proposed  platform  lends  itself  to  the 
induction  of  rocking  motions,  i.e.  the  generation  of  angular  motions  by 
translational  motions.  The  sketch  on  PI  13  explains  the  two  mechanisms 
through  which  these  rocking  motions  are  excited.  One  of  these  mech¬ 
anisms,  in  the  following  called  mechanism  C,  is  due  to  the  translational 
accelerat ion,  a,  operat ing  on  the  center  of  gravity  (CG)  of  the  primary 
reaction  mass  which  is  displaced  by  the  lever  arm  z  from  the  neutral 
plane.  In  this  way  a  torque  is  produced  which  leads  to  angular  motions 
about  CG.  The  other  mechanism,  called  mechanism  B  in  the  following,  is 
due  to  the  difference  between  the  translational  displacement  of  the 
ground  and  that  of  the  platform.  For  the  quantitative  analysis  of 
these  mechanisms  a  single  stage  instead  of  a  dual  stage  isolation 
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Lateral  ground  displacement  (in) 

Lateral  platform  displacement  (in) 

Lateral  platform  acceleration  (~A>*) 

2g  }  from  DD  7015-003 

2.8  in  from  p.E-3**  (Vol.  l) 

j  h  c pr 2 

Mass  of  Primary  mass;  kU.2  — ^ -  from  p.E-33 

lb  sec^ 

Mass  of  total  system;  68.7  — : -  from  p.  E-106 

in 

Moment  of  inertia  of  primary  mass;  3-25  *;  10**  lb  sec^  in 
Moment  of  inertia  of  total  system;  6.  A3  x  10*4  lb  sec^  in 


l=Em.d.  ;  m.  from  p.  E -  1  02 f f 

ill 

d.  from  p.  E 8 5 f f 

If  j  fc, 

Horizontal  spring  constant:  l.A2xl0  - — 

J  in 


Vertical  spring  constant:  3 -  1 3  x  10 


3  lb 


i  n 


calculated  from 
C - 55  and  0-5*4 
respect i vel y 


Rotational  spring  constant:  8.8x10°  lb  in  about  horizontal  axis:  S  = 
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system  was  assumed,  to  be  able  to  use  transfer  functions  from  the 
available  literature.  This  simplified  analysis  is  considered 
adequate  to  yield  order-of-magn i tude  information  on  the  effects.  The 
pertinent  coefficients  for  this  analysis  were  derived  from  the  trans- 
missibility  function  on  p.  C-5^  and  p.  C-55  under  the  assumption  that 
these  functions  were  generated  by  a  single  stage  system.  The  transfer 
function  for  mechanism  B  on  PI  14  is  transcribed  from  equation  (5)  of 
Ref.  5.  The  transfer  function  for  mechanism  C  is  straightforward.  The 
flow  chart  on  PI  15  shows  the  steps  to  be  taken  to  compare  the  various 
angular  motions  while  the  active  angular  stabilization  loop  is  not  in 
operation  but  the  horizontal  active  loop  is.  PI  16  is  a  compilation 
of  all  the  calculations  performed  and  PI  17  is  a  graphical  presentation 
of  the  results.  It  is  seen  that  the  various  motion  components  are 
essentially  of  the  same  order  of  magnitude.  There  is  no  evidence  of 
a  catastrophic  excitation  of  rocking  motions.  However,  the  analyst  of 
future  gyro  and  accelerometer  test  results  must  be  aware  that  such 
rocking  mechanisms  exist,  must  look  out  for  them  in  his  test  data,  and, 
if  they  are  found,  must  avoid  attributing  them  to  the  test  item. 

3.2  Instrumentat ion  Concept 

The  instrumentat ion  philosophy  of  subject  report  is  quite  inter¬ 
esting.  In  contrast  to  other  high  precision  stabilized  platforms  which 
use  two  different  types  of  instruments  to  be  able  to  measure  all  plat¬ 
form  motions  over  the  frequency  range  of  interest,  the  report  proposes 
to  use  only  one  type.  The  instruments  in  question  are  seismometers,  of 
which  4  will  be  provided  to  measure  motions  in  the  horizontal  and  4 
in  the  vertical  plane.  The  manufacturer  of  the  seismometers  is 
Sprengnether  Instruments,  Inc.,  and  the  type  numbers  are  S-5100-H  for 
the  horizontal  and  S-5100-V  for  the  vertical  instruments.  The  seis¬ 
mometers  are  of  the  spring-mass  type  and  the  deflection  of  the  proof 
mass  is  measured  either  by  a  coil  (velocity  pick-off)  or  a  capacitance 
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R  •  Horizontal  damping  force  coefficient:  3  •  1 7  x  lo'  -r^-  sec 
n  i  n  _  *  _  _ 

,  }  From  C-55 

Rv :  Vertical  damping  force  coefficient:  1.51x10'  yjy  sec  and 

H  :  Rotational  damping  force  coefficient  about  horizontal  axis: 

H,  =  R  =  4. 24x10^  lb  in  sec 
h  v 

Rocking  Mechanism  B 
0=  (Us-Ug)  (SZRh-Zkh) 

S2|  +  S(H,  +  Z2R,  )  +  (S,  +  Z2K,  ) 
t  n  h  h  n 

G  _  S(28  X  3- 17  X  101)  -  (28  x  1 .42  X  IQ4) _ _ 

AU  S2(6.43x  10A)  +  S(4.24  +  104+  7 -  84  x  10~2x  3 . 1  7  x  101)  +  (8.8  x  )06  x  7-84  x  102  x  1.42  x  10**) 

0  =  S(4.47  x  10~5)  -  (2  x  10'2) _ 

AU  S2  (3  •  23  x  10'3)  +  S  (3  -  38  x  10"3)  +  1 

Frequency  =2.8  Hz 
DC  Gain:  2x10  2 
Resonance  Peak  (relative):  16.8 


Rocking  Mechanism  C 
aM2  ‘  Z  =  0  [S2  1 2  +  SHh+  Sh] 

0  _  H2Z  .  0  _  (44.2  x  2.8)  _ 

a  S2|0  +  SH.  +  $  ’  r  S2  (3-24  +  10  )  +  S (4.24 +  10  )  +  (8.8  x  10  ) 
z.  n  n 

0 _  K4  x  10~5 _ 

a  Sz(3.68x  lo”3)  +  S  ( '• .  82  x  lo’3)  +  1 

Frequency =  2 . 62  Hz 

DC  Gain:  1.4x  10’5 

Resonance  Peak  (relative):  12.5 
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COMPARISON  OF  ANGULAR  PLATFORM  MOTION  COMPONENTS 
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bridge  (dl splacement  pick-off).  They  have  a  second  pick-off  coil  which 
can  be  used  to  adjust  the  damping  coefficient.  The  instruments  seem  to 
be  of  good  quality  and  are  surprisingly  inexpensive  so  that  the  high 
number  of  instruments  which  provides  a  certain  redundancy  is  justified. 

To  facilitate  the  understanding  of  how  the  seismometers  are 
proposed  to  be  used,  PI  18  gives  an  overview  of  their  pertinent  char¬ 
acteristics.  Their  transfer  function  is  of  second  order  and  the  expon¬ 
ent,  x,  of  the  numerator  function  depends  on  the  input  quantity  of 
interest  and  on  the  pick-off  from  which  the  output  is  taken.  With  x 
defined  it  can  be  shown  how  the  transfer  function  behaves  at  low 
frequencies  (below  the  natural  frequency)  and  at  high  frequencies 
(above  the  natural  frequency).  As  an  example,  if  linear  acceleration 
is  the  input  quantity  of  interest  and  the  output  is  taken  from  the 
velocity  pick-off  then  x  =  1  and  the  output  at  low  frequencies  is 
proportional  to  s  and  at  high  frequencies  to  l/s.  It  must  be  noted  that 
.the  instrument  is  capable  of  measuring  angular  displacement  because 
through  such  a  displacement  it  picks  up  a  component  of  the  gravity 
vector  and  the  resulting  measured  acceleration  can  be  interpreted  in 
■terms  of  angular  displacement. 

In  the  proposed  system  scheme  the  seismometers  measure  trans lat ional 
motions  in  their  respective  reference  planes  and  if  more  than  one  instru¬ 
ment  is  avaliable  in  each  direction  the  redundancy  is  used  to  establish 
■an  optimum  estimate.  Of  particular  interest  is  the  approach  taken  to 
measure  angular  platform  motions.  This  is  described  on  PI  19.  Using 

i 

as  an  example  the  tilt  motions  about  a  horizontal  axis,  the  sketch  of 
'  the  geometry  shows  that  one  horizontal  and  two  vertical  seismometers  are 
used  for  the  purpose.  At  low  frequencies,  the  angular  information  is 
taken  from  the  displacement  pick-off  of  the  horizontal  instrument,  which, 
according  to  the  table  in  PI  18  is  constant  in  this  region.  Of  course, 
the  instrument  is  also  sensitive  to  translational  accelerations  but 
these  are  negligible  in  this  frequency  range. 


SEISMOMETER  OVERVIEW 


General  Transfer  Function 


2  2 
S  +  2  ?  iii  S  +  w 
n  n 

For  Sprengnether  Seismometers  S-5100  H  and  S-5100  V 

=  0.003  x  6.28  =  1.89  xl0-2 
n  Sec 

5  is  adjustable  through  separate  pick-off  coil 


Input 

Velocity  pick-off 
K=KV 

Displacement  pick-off 

K=KP 

Linear 

acceleration 

x  =  1 

_ 

'  x  *  0 
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S 

mm 
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High 

1/S2 

Linear 

velocity 

X 

=2 

- 1 
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S2 
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S 

High 

1/S 

Linear 
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x  =  3 

x  —  2 
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S3 
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S 
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S2 
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Angular 

displacement 

_l 

x=  1 

x=  0 
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ANGULAR  MEASUREMENT 


Angular  information  output  a: 


At  low  frequencies: 

From  displacement  pick-off  of  H ^ :  otj  =  0*g:  constant  response 
At  high  frequencies: 

From  velocity  pick-off  of  V^and  V2  :  ap=vi~v2 

=  01  Ck1  +  k2) 

=  0s 1  (kj  +  k2) 

To  equalize  high  frequency  output,  a  filter  with  response  ^  ^  is  required 
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At  high  frequencies  the  angular  information  is  obtained  by 
differencing  the  velocity  outputs  of  the  two  vertical  seismometers. 

Each  instrument  measures  the  vertical  velocity  due  to  the  tilt  rate  as 
described  on  PI  19.  With  respect  to  tilt  rate  the  output  is  constant 
in  this  frequency  range,  but  if  the  tilt  angle  is  defined  as  the 
quantity  of  interest  then  the  output  is  proportional  to  s  and  a  filter 
with  an  intergrating  response  must  be  used  to  equalize  it.  Then,  if 
the  scale  factors  of  the  horizontal  seismometer  and  of  the  two  verti¬ 
cal  seismometers  are  matched  a  constant  t ransmi ss ibi 1 i ty  at  low  and 
high  frequencies  is  obtained.  There  are  some  deviations  in  the  trans¬ 
ition  zone  between  the  two  frequency  ranges  and  around  the  natural 
frequency  of  the  instrument.  This  is  discussed  in  detail  on  p.  D-15 
and  it  is  shown  that  these  deviations  are  tolerable  if  the  damping 
ratio  of  the  instruments  is  adjusted  to  0.1»-0.5. 

Now,  the  vertical  seismometers  measure  linear  vertical  velocity 
in  addition  to  the  velocity  induced  by  the  tilt  rate.  This  effect  is 
discussed  in  detail  on  PI  20.  In  theory,  this  undesired  velocity  out¬ 
put  is  eliminated  by  differencing  the  two  seismometers  outputs.  But  if 
the  scale  factors  of  the  two  instruments  are  not  exactly  equal  a  certain 
undesired  output  remains.  Based  on  the  quantitative  knowledge  of  the 
angular  and  linear  motions  of  the  platform,  the  last  column  of  the 
table  in  PI  20  shows  the  maximum  admissible  unbalance  of  the  two 
instruments.  The  requirements  are  high  but  tolerable  at  high  frequen¬ 
cies  but  increase  in  severity  with  decreasing  frequency.  Somewhat 
disturbing  is  the  fact  that  very  little  is  known  about  translational 
motions  in  the  frequency  range  below  0.2  Hz.  It  is  very  well  possible 
that  the  balance  requirements  become  even  more  severe  at  frequencies 
below  0.2  Hz  which  would  probably  make  the  proposed  angular  measurements 
scheme  impractical. 

It  is  therefore  absolutely  necessary  that  an  in-depth  analysis  of 
this  problem  area  is  performed  before  a  final  decision  with  respect 
to  instrumentation  is  made.  This  analysis  must  be  based  on  extensive 


Seismometer  Unbalance  Effect 


Actual  seismometer  outputs: 

V1  “  kl  (61 +  v) 

V2  -  k2  (-61  +  v)  velocity 

Vl~V2  =  01  (k1  +  k2) +v-(k1-k2) 

To  make  undesired  output  v(k^-k2)  equal  to  or  smaller  than  desired  output  01  (  kj+k2) 

Ak  01 
2k  =  v 

On  actively  stabilized  platform: 

0  =  0*T  ;  T  :  active  transmissibility  from  p.  D-30  0  Derived  from  PI  16 

P  1 

v  =  p;  a^_  taken  from  PI  15 

p  -  P 

r  ui 


1=66" 


f 

01 

in/Sec 

T 

0 

Op*  1 

in/ Sec 

1 

in/Sec^ 

- 

in/Sec 

V1 

VP 

30 

2.07  x 10"7 

3.2  x 10~3 

6.6  x  10"10 

4.6x10  8 

2.44 x  10~8 

2. 7  x  10~2 

20 

2.07  x  10-7 

1 .08  x  10-9 

8.4  x  10-7 

6.67 x  10~9 

1.6  x  10'1 

10 

2.07  x  10'7 

2.07  x  Iff 9 

8.6  x  10-7 

1.36  x  10~8 

1.52  x  10-1 

a 

2.07  x  10-7 

1.4  x  10~2 

2.9  x  10-9 

9.75  x  10"7 

. 

3.88  x  10~8 

7.47  x  10~2 

i 

2.07  x  10~7 

2  x  10'2 

4.1  x  10-9 

3.2  x  10“6 

5.1 x  10~7 

8  x  10~3 

0.2 

2.07  x  10~7 

2.8  x  10~2 

5.8  x  10-9 

3.1 x  10-6 

2.47 x  10“6 

2.35 x  10~3 
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in-situ  measurement  of  angular  and  linear  motions  and  on  a  critical 

evaluation  of  the  expected  stability  of  the  instruments.  The  fre- 

-1  -2 

quency  range  between  10  and  10  Hz  is  particularly  problematic 
because  it  involves  ground  motions  which  are  caused  by  the  so-called 
microseisms  which  in  turn  are  due  to  ocean  waves  pounding  on  the 
continental  shelf  and  which  are  not  fully  understood.  The  discus¬ 
sion  on  page  3  of  Appendix  Bwill  lead  to  a  better  understanding 
of  this  phenomenon. 

3.3  Active  Loop  Analysis 

While  introducing  some  sound  and  interesting  ideas  such  as 
using  only  one  type  of  instruments  and  using  a  digital  computer  in 
the  control  loop,  the  presentation  of  the  Design  Analysis  for  the 
Active  Control  System  in  chapter  D  of  subject  report  is  rather  sloppy. 
This  is  so  because  in  many  cases  unit  dimension  information  is  missing 
or  erroneous,  symbols  are  not  defined  or  confusing,  indices  are  mixed 
up,  equations  are  not  numbered  and  in  some  cases  incomplete,  and  the 
validity  of  some  results  is  suspect.  The  major  deficiency  is  the 
lack  of  emphasis  given  to  the  analysis  of  rotational  (angular)  motions 
which  is  particularly  deplorable  because  these  motions  are  the  main 
contributor  to  test  errors  as  demonstrated  in  par.  2  above.  The 
analysis  conducted  by  this  author  therefore  will  concentrate  on 
angular  motion  to  fill  this  gap. 

The  first  task  is  to  calculate  the  passive  angular  transmissi- 
bility  function  which  can  not  be  found  anywhere  in  the  report.  PI  21 
presents  the  pertinent  transfer  function  which  is  structured  in  ac¬ 
cordance  with  the  translational  t ransmi ss i b i 1 i ty  function  on  p.  D~3. 
This  is  done  despite  the  fact  that  it  does  not  agree  with  the 
equations  of  rotational  motion  on  p.  D-5.  But  these  are  suspect 
because  they  do  not  agree  with  the  block  diagram  on  p.  D-7  which  seems 
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PASSIVE  ANGULAR  TRANSFER  FUNCTION 


V 


V 


S2  •  E(3)  4-  S  ■  E(2)  +  F(n _ 

S4  •  F(5)  +  S3  •  F (4)  +  S2  •  F(3)  +  S  •  F(2)  +  F(l) 
Angular  displacement  of  primary  reaction  mass 

Angular  displacement  of  ground 


E(3)=  Bx  •  B2 
E(2)=  K1B2  +  K2B1 
F(l)=  KjKg 


F(2)=  K1B2  +  K2B1 

F(3)“  I2Ki  +  12K2  +  BiB2  +  IiK2 
F(4)-  I^  +  IjBj+Vz 
F(5)-  Ijl2 


*2  :  Moment  of  inertia  of  primary  reaction  mass:  8.13xl03  [lb  in  Sec2] 

Ij  :  Moment  of  inertia  of  secondary  reaction  mass:  7.43xl03  [lb  in  Sec2] 

"2 
K 


Spring  constant  of  upper  springs:  3.42xl06  [lb  in] 
]  :  Spring  constant  of  lower  springs:  5.12  x  106  [lb  in] 


b2:  Torque  coefficient  of  upper  springs: 
B, :  Torque  coefficient  of  lower  springs: 


}  9.7  x  10  [lb  in  Sec] 


NOTE:  All  coefficients  are  defined  per  spring  and  not  for  the 
whole  system. 
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to  represent  the  physical  reality  quite  well.  PI  21  further  gives 
the  definition  of  the  coefficients  of  the  transfer  function  for 
tilt  motions  about  a  horizontal  axis  and  a  quantitative  list  of  the 
pertinent  physical  coefficients.  They  are  derived  from  the  coeffi¬ 
cients  which  lead  to  the  translational  transmiss ib i 1 i ty  function  on 

2 

p.  C- 54  by  multiplying  the  spring  related  constants  by  d  ,  where  d 
is  the  distance  of  the  spring  from  the  center  line  of  the  platform. 

A  list  of  rotational  coefficients  given  in  Table  0-1  on  p.  0-19  is 
highly  suspect  and  was  not  used  because,  besides  erroneous  unit 
dimension  information,  it  gives  different  values  for  the  inertia  in 
the  two  horizontal  directions  which  can  not  be  explained  because  the 
platform  structure  is  symmetrical. 

Furthermore,  it  looks  like  some  of  the  coefficients  are  defined 
for  the  whole  system  while  others  are  on  a  "per-spr 1 ng"  basis.  The 
top  curve  on  PI  22-  shows  the  passive  t ransmi ss i bi 1 i ty  function  for 
tilt  motions  about  a  horizontal  axis.  It  will  be  noted  that  it  is 
quite  similar  to  the  vertical  translational  transfer  function  of 
p.  C- 54,  as  must  be  expected.  In  particular  it  shows  the  steep  4th- 
order  drop-off  at  high  frequencies,  a  factor  which  is  important  for 
the  following  discussions. 

For  the  analysis  of  the  active  angular  control  system,  PI  23 
shows  first  the  transfer  function  describing  the  response  of  the 
platform  angle  to  torques  applied  to  the  platform  by  the  electro¬ 
magnetic  actuators.  The  equation  is  formulated  in  accordance  with 
the  equation  for  translational  motions  on  p.  D-5.  The  pertinent 
additional  coefficients  are  defined. 

The  dependence  of  the  stabilization  torque  on  the  platform 
angle,  the  missing  link  for  loop  closure,  is  described  by  the  next 
equation.  Its  simple  form  is  justified  for  the  following  reason: 

As  explained  in  PI  13  and  PI  19,  the  angular  transmi ss i b i 1 i ty  of  the 
seismometers  with  equalizing  filter  in  the  critical  high  frequency 
range  is  constant.  The  electromagnetic  actuators  are  expected  to 
have  a  bandwidth  of  20  Hz  (see  p.  24),  so  they  do  not  contribute 
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ACTIVE  ANGULAR  CONTROL  LOOP 

?2  m  S2  •  Il  +  S  (Bj  +  B2)  +  (K:  +  K2) 

S4  •  F(5)  +  S3  •  F(4)  +  S2  •  F(3)  +  S  •  F (2)  +F(1) 

T:  Torque  applied  to  primary  reaction  mass  by  electromagnetic  actuators.  [lb  in] 

-  =  K„  +  SKi 
2  0  0 

K0:  Position  feedback  coefficient  [lb  in] 

K^:  Rate  feedback  coefficient  [lb  in  sec] 


CLOSED  LOOP  TRANSFER  FUNCTION 

0D  [S2  •  E  (3)  +  S  •  E(2)  +  F  (1)  ]  -  02  (KQ+SK^)  [S2I1  +  S  (BJ  +  B2)  +  (Kj  +  1^)  ] 

2  S4  •  F(5)  +  S3  •  F (A j  +  S2  •  F(3)  +  S  •  F(2>  +  F(l) 

@2  2 

— :  Numerator  function:  S  •  E(3)  +  S  •  E(2)+F(l) 

Denominator  function:  S  *  F(5)  +  S2  *  (F(4) +  }  +  S2  •  (F(3)  +  Kg  I ^  (B^  4-  B2) } 

+  s  ■  {f(2)  +  k0(b1  +  b2)  +  k^  (k1  +  k2))+{f(1)  +  k0  (Kj  +  Ko) ) 

0^:  Angular  Ground  Motion 
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any  dynamic  terms  in  the  frequency  range  of  interest  either.  Thus 
it  is  possible  to  describe  the  feedback  as  the  sum  of  position  and 
rate  feedback  gains.  An  acceleration  feedback  was  not  considered  in 
this  analysis  nor  were  there  any  additional  compensation  terms  intro¬ 
duced  . 

Finally,  PI  23  presents  the  complete  closed  loop  transfer  func¬ 
tion.  Thanks  to  the  simplified  form  of  the  feedback  term,  the 
closed  loop  equation  is  manageable  for  evaluation  by  pocket  calcula¬ 
tor  and  two  results  are  shown  on  PI  22. 

These  results  deserve  thorough  discussion.  Concentrating  first 
or  the  middle  one  it  is  seen  that  it  differs  significantly  from  the 
active  angular  transmi ss i b i 1 i ty  of  subject  report  on  p.  D-30  which 
was  the  basis  for  all  analyses  in  previous  chapters  of  this  evalu¬ 
ation.  While  the  two  agree  well  at  low  frequencies  (the  feedback 
coefficients  were  selected  accordingly)  there  are  wide  discrepancies 
at  high  frequencies.  In  particular,  the  steep  Ath  -  order  drop-off 
is  missing  at  the  MAC  function.  This  is  important  because  the  higher 
attenuation  at  high  frequencies  which  contribute  significantly  to  the 
test  errors  indicates  that  the  active  angular  stabilization  system 
actually  is  better  than  anticipated  from  the  results  of  subject 
report.  Of  course,  the  question  arises:  Which  of  the  two  functions 
is  valid?  The  fact  that  the  high  frequency  active  response  coincides 
well  with  the  passive  response  for  translational  motions  (see  PI  2) 
indicates  that  the  same  must  be  expected  for  angular  motions.  It 
does  not  make  any  sense  if  the  active  response  is  inferior  to  the 
passive  response  at  these  frequencies,  because  then  the  active  sys¬ 
tem  instead  of  stabilizing  would  be  destabilizing.  There  is  no 
evidence  that  subject  report  uses  any  additional  compensation  terms. 
The  only  difference  between  the  two  approaches  is  with  respect  to 
the  feedback  terms.  The  formulas  on  top  of  p.  D- 1 8  show  the  feedback 
term  used  by  MAC.  In  addition  to  position  and  rate  feedback  it  uses 
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acceleration  feedback  and  a  low-pass  filter  with  31.8  Hz  natural 
frequency  (see  Table  D-2  on  p.  D - 1 9 ) •  The  term  is  somewhat  suspect 
because  it  goes  to  0  i f  there  is  no  acceleration  feedback.  In  any 
case,  use  of  this  feedback  term  can  not  explain  the  slow  drop-off  at 
high  frequencies.  It  must  then  be  concluded  that  our  transfer  func¬ 
tion,  though  based  on  certain  simplifications,  is  valid  while  the 
MAC  function  is  not,  because  its  character  can  not  be  explained 
with  the  information  available  at  this  point. 

It  is  of  importance  to  point  out  that  the  loop  is  stable  under 
the  feedback  conditions  selected.  This  suggests  to  try  an  even 
tighter  feedback.  The  lower  curve  on  PI  2  shows  one  result.  Its 
general  attenuation  level  is  about  one  order  of  magnitude  higher 
than  before  and  the  loop  is  still  stable.  However,  there  is  indica¬ 
tion  that  the  stability  limit  is  being  approached  because  the 
resonance  peak  at  5.2  Hz  is  more  pronounped.  There  was  no  time  to 
determine  the  stability  margins  by  open  loop  analysis. 

In  conclusion,  it  can  be  stated  that  there  is  strong  evidence 
that  the  active  angular  control  system  is  better  than  expected  by 
the  contractors  analysis.  This  is  due  to  the  fact  that  there  is 
significantly  higher  attenuation  at  high  frequencies  and  that  the 
feedback  loop  can  be  tightened  without  endangering  stability.  The 
potential  improvement  is  estimated  at  one  to  one  and  a  half  order 
of  magnitude  which  would  make  it  possible  to  reach  the  test  goals 
spelled  out  on  PI  8  and  PI  10.  However,  caution  is  indicated 
because  simplifications  were  used  in  the  present  analysis  and  an 
in-depth  computer  analysis  is  required. 

.  Conclusions  and  Pecommendat  ions  : 

1.  The  passive  isolation  system  is  basically  sound  and  in  combin¬ 
ation  with  the  active  system,  as  described  in  subject  report,  is 


compatible  with  the  state  of  the  art.  However,  there  Is  strong 
evidence  that  the  active  stabilization  loop  actually  is  better  and 
can  be  even  more  improved.  In  this  way  it  seems  possible  to  even¬ 
tually  provide  a  capability  to  test  gyros  with  a  drift  rate  of 

-6  "8 
10  (deg/hr)rms  and  accelerometers  with  a  noise  level  of  10  grms. 

2.  There  is  a  certain  crosscoupi ing  between  translational  and 
rotational  motions  due  to  the  excitation  of  rocking  motions.  At 
the  given  disturbance  input  levels,  these  rocking  motions  are  of 
the  same  order  of  magnitude  as  the  natural  rotational  motions  and 
as  such  are  not  dangerous.  However,  the  analysis  performed  used 
certain  simplifications  and  it  is  necessary  to  perform  an  in-depth 
computer  analysis  when  the  final  system  constants  are  known  and 
more  information  on  the  disturbance  inputs  at  the  test  site  will 
be  avai lable. 

3.  A  similar  crosscoupling  between  translational  and  rotational 

motions  is  due  to  the  gain  unbalance  of  the  two  seismometers  which 

measure  angular  motions.  The  present  'analys i s  indicates  that  the 

relative  difference  of  the  gain  factors  must  not  exceed  a  few  parts 
-k 

in  10  to  avoid  that  the  induced  angular  motion  is  larger  than  the 
natural  angular  motion.  This  is  already  a  requirement  on  the  bor¬ 
der  line  of  realizability.  It  is  therefore  necessary  that  an  in- 
depth  analysis  be  performed,  based  on  the  most  reliable  measure¬ 
ments  of  the  input  disturbances  at  the  test  site  and  a  critical 
review  of  instrument  stability,  before  a  final  decision  is  made  on 
the  scheme  of  using  two  balanced  seismometers  for  measuring 


k.  From  the  foregoing  analysis  it  is  evident  how  important  the 

knowledge  of  actual  angular  and  translational  motions  at  the  test 

site  is.  In  particular,  it  is  necessary  to  provide  more  information 
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at  least  down  to  frequencies  of  about  10  Hz  to  cover  the  range  of 
the  so-called  microseisms.  The  interaction  of  translational  and 
rotational  motions  in  this  range  depend  heavily  on  the  test  site. 

It  is  therefore  necessary  to  perform  an  extended  series  of  measure¬ 
ments  at  the  test  site,  using  instrumentation  which  clearly  distin¬ 
guishes  between  rotational  and  translational  motions. 

5-  While  not  directly  a  part  of  this  evaluation,  it  is  suggested 
that  the  analysts  of  future  tests  performed  on  the  platform  develop 
a  deep  understanding  of  the  platform  characteristics.  This  will 
enable  them  to  recognize  events  which  are  caused  by  the  platform 
and  not  by  the  test  item  and  to  eliminate  them  from  test  data. 
Examples  are  rocking  motions  of  the  platform.  Related  to  this 
are  events  like  airplane  overflights  or  sled  runs  the  time  of 
occurence  of  which  can  easily  be  identified  and  their  effects  on 
the  test  data,  if  any,  can  be  eliminated. 

6.  A  further  improvement  of  test  accuracy  can  be  achieved  if  the 
test  analyst  critically  determines  what  frequency  ranges  are  of 
actual  interest  to  the  item  under  test.  As  an  example,  a  gyro 
under  test  may  not  have  a  quick  reaction  requirement  and  conse¬ 
quently  corresponding  higher  frequencies  can  be  filtered  out  from 
the  test  data.  This  can  lead  to  a  substantial  reduction  of  noise 
1 evel . 


Appendix  A. 

Discussion  of  Data  Presentation  Format 

The  variables  we.  have  to  deal  with,  rotational  and  transla¬ 
tional  ground  motions  and  their  derived  quantities,  are  random 


processes.  The  frequency  components  of  such  processes  do  not  have 
defined  phase  relationship  and  thus  the  only  possible  way  of  their 
spectral  representation  is  by  the  Power  Spectral  Density  Functions 
(PSD)  or  Power  Spectra, which  have  a  dimension  of  sec^/Hz,  using  a 
rotational  quantity  as  an  example.  For  a  detailed  discussion  of 
this  see  Ref.  6.  These  power  spectra  make  it  possible  to  determine 
the  over-all  rms  value  of  the  variable  in  a  given  frequency  band  by 
integrating  the  power  spectrum  over  this  frequency  band.  In  this 
way  the  total  "power"  in  the  band  is  obtained  and  by  taking  the 
square  root  the  rms  value.  In  the  Second  Design  Review  Meeting  of 
8  Mar  80,  it  was  agreed  to  use  a  "rms  spectrum"  or  "amplitude  spec¬ 
trum"  for  the  presentation  of  the  data,  in  our  example  having  a 
dimension  of  Secrms/Hz.  This  can  be  obtained  by  taking  the  square 
root  of  the  power  spectrum  at  each  frequency.  This  is  perfectly 
legal  and  physically  meaningful.  Instead  of  saying:  my  variable  has 
so  much  "power"  in  a  frequency  interval  of  1  Hz  it  is  stated  now 
that  it  has  such  and  such  rms  value  in  this  interval  of  1  Hz. 

However,  there  is  a  possible  trap:  It  is  not  possible  to  determine 
the  over-all  rms  value  of  a  variable  by  integrating  the  rms  spectrum 
over  frequency.  To  determine  the  overall  rms  value  it  is  necessary 
first  to  obtain  the  power  spectrum  by  squaring  the  rms  spectrum  point 
for  point,  then  integrating  the  power  spectrum  over  frequency  to 
obtain  total  power  and  finally  by  taking  the  square  root  obtainning 
the  over-all  rms  value. 

Appendix  B 

See  attached  paper  titled,  A  Proposed  Model  for  the  Rotation¬ 
al  Motions  of  a  Test  Pad. 
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Abstract 

Based  on  pad  motion  measurements  at  various 
test  sites,  the  model  is  formulated  as  the  power 
spectral  density  function  of  angular  deflections. 
Basically,  it  consists  of  a  random  walk  process 
which  is  well  defined  over  10  decades  of  frequency. 
Superimposed  are  two  peak  areas,  one  of  them  re¬ 
flecting  diurnal  fluctuations  and  the  other  motions 
due  to  microseisms.  The  properties  of  these 
processes  and  their  estimated  tolerances  are  dis¬ 
cussed.  However,  the  mechanism  of  the  fluctuations 
due  to  microseisms  is  not  yet  fully  understood  and 
additional  experimental  and  theoretical  work  is 
solicited.  Derived  from  the  model  are  the  expected 
uncertainty  in  the  orientation  of  an  azimuth 
reference  monument  as  a  function  of  time  after 
calibration  and  the  angular  rates  which  will  be 
sensed  by  a  gyro  under  test. 

I,  Introduction 

The  ever  increasing  accuracy  of  gyroscopes 
makes  the  testing  of  these  instruments  without 
undue  interference  by  external  perturbations 
increasingly  difficult.  One  source  of  perturba¬ 
tions  which  is  of  particular  significance  is  the 
rotational  motion  of  the  pad  on  which  the  gyro  is 
tested,  a  motion  which  is  represented  by  its  three 
components:  North-South  (NS)  tilts,  East-West  (EW) 
tilts,  and  changes  in  azimuth  orientation.  The 
concern  with  this  motion  environment  is  certainly 
not  new.  Extensive  measurements  of  the  motion 
components  have  been  made  by  many  test  agencies  and 
various  approaches  to  obtain  control  over  these 
motions  are  under  way.  However,  what  is  missing  at 
this  time  is  a  quantitative  understanding  of  the 
functional  character  of  this  motion  environment. 
Availability  of  such  an  understanding,  in  other 
words  the  existence  of  a  "model"  would  facilitate 
the  optimization  of  test  procedures  in  the  follow¬ 
ing  areas:  Calibration  of  the  azimuth  orientation 
of  a  test  pad,  compensation  of  tilt  and  azipiuth 
motions  and  optimum  filtering  of  the  gyro  outputs 
during  a  test. 

It  is  therefore,  the  purpose  of  t^is  paper  to 
review  and  compile  data  which  were  taken  at 
representative  test  sites  with  the  goal  of  compar¬ 
ing  magnitudes  and  detecting  common  trends.  As  a 
result,  a  quantitative  model  will  be  formulated  and 
its  uncertainties  will  be  discussed.  Thus  the 
modeling  process  is  purely  empirical  and  inductive, 
concentrating  on  the  phenomena  which  are  of  practi¬ 
cal  interest  to  the  test  engineer  and  attempts  of 
a  physical  explanation  of  the  model  will  be  very 
limited.  However,  the  surprising  agreement  with  a 
theoretical  model  which  was  derived  deductively 
many  years  ago  will  be  noted  and  the  ensuing  possi¬ 
bilities  for  a  physical  explanation  of  the  observed 
phenomena  will  be  discussed.  Finally,  basic  con¬ 
siderations  concerning  the  application  of  the  pro¬ 
posed  model  will  be  presented. 

*  Associate  Fellow  AIAA 


The  ultimate  purpose  of  the  paper  is  an  appeal 
to  the  Geokinetic  Community  to  create  a  data  pool 
which  is  much  more  comprehensive  than  the  one 
which  is  presently  available.  This  then  will  allow 
formulating  a  model  with  much  higher  confidence 
than  is  presently  possible. 

II.  Compilation  of  Data  and 
Formulation  of  the  Model 

An  overview  of  available  rotational  motion 
data  indicates  that  deterministic  trends  are  either 
non-existent  or  very  weak.  The  analytical  tool  for 
describing  such  data  thus  has  to  be  of  statistical 
nature  and  for  the  purposes  of  this  paper  the  tools 
of  time-series  analysis  will  be  used.  In  particular 
power  spectral  density  functions,  in  short,  power 
spectra,  will  be  used  to  describe  the  measured  data. 
This  then  allows  using  the  related  descriptors, 
namely,  cross-correlation  spectra,  transfer  functions 
and  coherence  functions  as  powerful  tools  for 
further  investigations.  Power  spectral  analysis  is 
not  new  but  has  come  into  widespread  use  only  fairly 
recently  after  the  full  development  of  the  Fast 
Fourier  Transform  and  the  availability  of  mini¬ 
computers  made  the  on-line  and  near  real-time  com¬ 
putation  of  power  spectra  a  practical  possibility. 
Reference  1  describes  this  analytical  method  in 
detail  and  discusses  its  limitations. 

The  frequency  range  of  a  power  spectrum 
describing  the  angular  motions  which  are  of  interest 
for  gyro  testing  is  extremely  wide,  extending  from 
about  100  Hz  to  10~8  Hz,  a  span  of  10  orders  of 
magnitude.  The  upper  limit  is  due  to  the  necessity 
of  determining  gyro  noise  within  the  pass-band  of  a 
platform  stabilization  loop,  the  lower  one  is 
connected  with  the  interest  in  the  behavior  of 
azimuth  reference  monuments  and  in  the  long-term 
stability  measurements  of  gyros.  As  will  be  shown, 
it  is  possible  to  formulate  a  model  in  terms  of  an 
angular  deflection  power  spectrum  which  is  meaning¬ 
ful  over  this  wide  frequency  range.  The  quantity 
of  interest J  angular  deflection  power  density, 
covers  a  span  of  20  orders  of  magnitude.  Quite 
certainly  a  phenomenon  extending  over  such  an 
extremely  wide  range  of  magnitudes  cannot  be  ob¬ 
served  by  a  single  type  of  instrument  or  one  single 
method.  It  is  rather  necessary  to  perform  measure¬ 
ments  within  restricted  frequency  bands  using  in¬ 
struments  and  methods  which  are  suitable  for  the 
respective  range. 

Representative  results  of  measurements  made 
in  various  frequency  ranges  are  presented  in  the 
following  figures.  The  data  used  are  extracted 
and  transcribed  from  the  sources  which  are  identi¬ 
fied  in  the  figures.  To  make  the  conyilation  of 
many  data  in  one  figure  possible  and  still  avoid 
overcrowding,  the  connecting  lines  between  data 
points  are  omitted.  Furthermore,  only  such  data 
points  were  transferred  from  the  original  source 
which  are  required  to  reflect  realistically  the 
trends  and  uncertainties  of  a  given  spectrum. 
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Figure  1  covers  the  two  decades  from  1  Hz 
to  100  Hz.  The  measurements  are  made  with  arrays 
of  seismometers.  In  this  way  it  is  possible  to 
measure  all  three  components  of  angular  motion,  NS 
tilt,  EW  tilt,  and  fluctuations  in  azimuth.  An 
important  phenomenon  in  this  frequency  range  are 
spikes  which  are  clearly  due  to  cultural  activities, 
such  as  the  spikes  at  60  Hz  and  near  30  Hz.  Figure 
2  shows  the  two  tilt  components  in  the  range  from 
10*4  Hz  to  1  Hz.  The  measurements  are  made  with 
tilt  meters,  azimuth  information  is  not  available 
in  this  range.  Of  particular  interest  is  a  peak 
centered  at  about  0.15  Hz.  It  is  due  to  micro- 
seismic  activities  and  will  be  discussed  in  detail 
later.  Figure  3  covers  the  range  from  about  10~6 
Hz  to  10*2  Hz.  Again,  the  measurements  are  made 
with  tiltmeters  and  azimuth  information  is  not 
available. 

Finally,  Figure  4  presents  azimuth  data  in  the 
range  from  about  10"8  Hz  to  10-6  Hz.  Tilt  informa¬ 
tion  in  this  range  is  not  available.  The  data  are 
derived  from  the  calibration  of  the  azimuth  orienta¬ 
tion  of  reference  monuments  with  respect  to  Polaris. 
This  type  of  data  is  usually  not  very  well  suited 
for  power  spectral  analysis  because  the  observa¬ 
tions  are  normally  unequally  spaced  in  time.  How¬ 
ever,  the  spectrum  shown  was  derived  from  a  set  of 
measurements  which  were  fairly  regularly  spaced  at 
intervals  of  two  months,  and  thus  it  can  be  consid¬ 
ered  meaningful.  But  even  if  only  irregularly 
spaced  data  are  available,  it  is  still  possible  to 
make  a  reasonable  low-resolut ion  estimate  of  power 
density.  Assume  that  n  measurements  are  available, 
irregularly  spaced  over  an  observation  interval  T. 
Then  the  variance  of  the  measurements  can  be  calcu¬ 
lated  after  extraction  of  the  mean  and  linear  trend 
from  the  data.  This  variance  is  the  total  power  in 
the  frequency  band  between  1/T  Hz  and  n/2T  Hz  and  a 
mean  power  density  in  the  theoretical  frequency  band 
of  width  1/T  (n/2  -  1)  Hz  can  be  calculated.  Some 
results  of  this  type  are  shown  in  Figure  4.  This 
method  admittedly  is  crude  and  has  its  shortcomings, 
in  particular  with  respect  to  the  definition  of  the 
upper  frequency  limit.  However,  it  i,  considered 
adequate  for  an  order- of- magn i tude  check  of  spectral 
density  of  a  given  series  of  irregularly  spaced 
observations. 

To  obtain  a  comprehensive  overview  of  the 
character  of  the  rotational  motion  under  study,  the 
information  from  Figures  1  -  4  is  compiled  in 
Figure  5.  Again,  only  such  data  are  transferred 
which  are  necessary  and  sufficient  to  provide  ade¬ 
quate  trend  and  dispersion  information.  A  first 
inspection  of  the  figure  shows  that  the  individual 
spectra  follow  quite  closely  a  l/f2  law,  represented 
by  line  A,  over  the  wide  frequency  range  from  10'8 
Hz  to  102  Hz.  Considering  that  the  data  used  are 
derived  from  three  different  motion  components  and 
from  measurements  made  at  different  locations  at 
different  times,  the  dispersion  about  this  line  is 
quite  reasonable.  This  line  describes  the  under¬ 
lying  statistical  process  of  the  proposed  model  and 
is  defined  by  the  following  equation. 
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The  uncertainty  of  ±1  order  of  magnitude  may  seem 
large  but  is  considered  reasonable  for  a  phenomenon 
whose  magnitude  range  of  interest  extends  over  20 
decades. 


The  underlying  statistical  process  then  i 
random  walk.  There  is  some  evidence  that  this 
process  is  actually  an  "exponentially  correlated 
process",  the  spectral  density  of  which  would 
remain  constant  below  a  certain  breakpoint  fre¬ 
quency.  But  the  available  experimental  data  are 
too  scarce  to  make  a  reliable  quantitative  estimate 
of  such  a  breakpoint  frequency,  or,  in  other  words, 
of  the  correlation  time  of  the  process.  The  data 
compiled  in  Figure  5  indicate  that  for  observation 
periods  up  to  one  year,  or  frequencies  down  to  10-8 
Hz,  the  assumption  of  a  random  walk  is  a  valid 
representation  of  the  underlying  statistical  process 
of  the  proposed  model.  It  should  be  noted  that  the 
assumption  of  a  random  walk  character  for  the 
angular  displacements  leads  to  a  "white"  character 
for  the  angular  rates,  a  result  which  is  intuitively 
appealing. 

The  derivation  of  this  basic  model  implies  that 
within  the  given  tolerances,  its  three  components 
are  equal  and  that  it  is  independent  of  the  time  and 
place  of  the  observation.  A  review  of  existing 
literature  reveals  strong  indications  that  these 
assumptions  are  reasonable.  Line  B  in  Figure  5 
describes  a  model  predicted  by  Weinstock(2)  under 
the  assumption  that  the  rotational  phenomenon  is  due 
to  a  vertical  displacement  wave  propagating  in  an 
elastic  half-space  with  a  propagation  velocity  of 
about  two  miles/second.  In  accordance  with  the 
reference  the  rotational  displacement  spectrum  is 
then  given  as: 
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In  the  frequency  range  from  10  Hz  to  10  Hz 
Weinstock  assumes  constant  "white"  vertical  accel¬ 
eration  power  density  and  his  prediction  agrees 
surprisingly  well  with  the  proposed  model  not  only 
trendwise  but  also  magnitudewise.  That  the  assump¬ 
tion  of  a  white  character  of  the  translational 
acceleration  spectrum  is  reasonable  is  confirmed  by 
quite  recent  measurements  published  in  Reference  3. 
The  portions  of  these  measurements  pertinent  to 
the  present  discussion  cover  the  frequency  range 
from  6.6  x  10"^  Hz  to  3.3  x  10‘2  Hz.  In  addition  to 
confirming  the  character  of  the  spectrum  in  this 
frequency  range,  the  measurements  indicate  that  the 
spectral  magnitudes  of  the  three  motion  components 
are  equal  to  within  ±20%.  (Figure  9  of  Reference 
3),  that  the  observed  magnitudes  are  equal  within 
125%  independent  of  the  time  of  the  observation 
(Figure  10  of  Reference  3)  and  that  they  are  equal 
within  ±20%  independent  of  the  place  of  observa¬ 
tion.  (Figure  12  of  Reference  3). 


Superimposed  to  the  underlying  random  walk 
mechanism  which  is  described  above  are  several  peak 
areas.  The  one  which  turns  out  to  be  most  important 
is  the  one  which  centers  at  the  frequency  of  1.15  x 
10-5  Hz,  the  period  of  which  is  one  day.  The  strong 
influence  of  the  diurnal  cycle  is  well  known  and 
commonly  observed.  However,  it  is  also  known  that 
the  resulting  angular  motions  are  not  deterministic 
and  not  clearly  correlated  to  external  influences, 
such  as  temperature.  Apparently,  the  combination 
of  physical  phenomena  which  affect  the  angular 
motions  of  the  diurnal  cycle  do  not  fluctuate  peri¬ 
odically  but  in  a  random  fashion.  The  modeling 
concept  which  seems  to  explain  the  observations 
quite  realistically  is  that  of  a  sinusoidal  carrier 
frequency,  the  basic  diurnal  rhythm,  which  is  am¬ 
plitude  modulated  by  random  events.  The  spectrum 
resulting  from  such  a  process  can  be  described  by 
the  following  equation: 
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:  Spectral  density  of  resulting  angular 
displacements. 

♦D  :  Spectral  density  at  1.15  x  10’  Hz,  the 

o  carrier  frequency. 

fD  *  1.15  x  10' 5  Hz 

B  Bandwidth  of  the  modulating  random 

fluctuations  (Hz). 

n  Empirical  exponent. 

A  review  of  the  available  data  indicates  that  they 
can  be  represented  most  adequately  if  one  assumes 
B  «  fp  and  n  =  4.  Physically  this  means  that  the 
spectrum  of  the  modulating  random  forces  contains 
primarily  frequencies  with  periods  longer  than  one 
day  and  falls  off  sharply  at  higher  frequencies. 

The  resulting  spectrum  shown  in  Figure  5  is  calcu¬ 
lated  under  the  assumptions  that  is  100  times 

o 

larger  than  the  spectral  density  of  the  random  walk 
at  this  frequency,  that  f^/B  =  1,  and  that  r.  -  4. 

Another  significant  peak  in  the  composite 
spectrum  is  centered  close  to  10'1  Hz  and  is  known 
to  be  due  to  microseisms  generated  by  ocean  waves. 

It  is  observed  all  over  the  continent  if  the  angu¬ 
lar  measurements  are  made  with  tiltmeters.  However, 
it  has  not  been  observed  with  gyros  or  with  other 
instruments  which  sense  angular  motions  only  and 
are  not  sensitive  to  translational  accelerations 
like  tiltmeters.  An  example  of  a  nongvroscopic 
instrument  of  this  type  is  the  fused  silica  angular 
motion  sensor  described  by  S.  Okubo  in  Reference  4. 
This  phenomenon  of  apparent  nonobservability  of  the 
microseisms  with  instruments  which  are  sensitive  to 
rotational  motions  only  seems  to  contradict  the 
expectation  of  coupling  between  translational  and 
rotational  motions  which  was  described  in  equation 
(2).  However,  Reference  3  points  out  that  the 
earth  noise  within  the  microseismic  peak  is  of 
oceanic  origin  but  in  the  adjoirir.g  V  .-or  frequency- 
range,  the  range  which  the  proposed  model  describes 
as  a  random  walk,  it  is  of  atmospheric  origin.  That 
the  propagation  characteristics  of  the  noise  differ 
substantially  between  the  two  ranges  is  shown  in 
curve  A  of  Figure  6.  This  curve  is  transcribed  from 
Figure  4  of  Reference  3.  j«  show:  the  squared  coher¬ 
ency  function  of  two  simultaneous  vertical  transla¬ 
tional  displacement  measurements  made  with  !wc  seis¬ 
mometers  placed  at  a  dist  ince  of  :3  -  m  rsr p.  ..x  mate!) 


543  m  belpw  the  local  surface.  While  the  coherence 
is  quite  strong  in  the  microseismic  peak  range,  it 
falls  off  sharply  and  significantly  toward  lower 
frequencies.  Another  indication  of  the  difference 
of  the  microseismic  mechanism  from  the  one  in  adjoin¬ 
ing  frequency  ranges  is  given  by  curve  B  in  Figure  6. 
It  is  transcribed  from  Figure  3  of  Reference  10  and 
shows  an  example  of  power  attenuation  of  the  noise 
with  depth  below  the  surface.  While  the  attenuation 
within  the  microseismic  peak  range  is  small,  it 
increases  significantly  at  the  transition  to  higher 
frequencies.  Thus,  with  the  experimental  evidence 
presently  available,  it  must  be  concluded  that  the 
spectral  peak  about  0.1  Hz  is  not  truly  a  component 
of  the  rotational  motion  environment  and  thus  it  will 
not  be  included  in  the  presently  proposed  model. 

Other  spectral  peaks  due  to  natural  causes 
which  one  has  to  expect  are  located  at  the  frequency 
the  period  of  which  is  one  year,  and  another  one 
with  a  period  of  a  half  day.  Due  to  the  extremely 
long  periods  of  observation  required  to  describe  the 
annual  peak,  very  little  experimental  material  is 
presently  available.  There  is  some  evidence  of  such 
a  peak  in  Figure  4,  but  the  uncertainties  of  this 
one  spectrum  are  too  high  to  allow  a  serious  model¬ 
ing  effort.  The  half-day  peak  is  occasionally 
observable  but  is  normally  submerged  in  the  envelope 
of  the  diurnal  peak. 

At  the  high-frequency  end  of  the  spectrum  very 
high  and  narrow  spectral  peaks  can  be  observed.  (See 
Figure  1)  As  stated  before,  they  are  due  to  human 
activities  and  can  be  easily  identified  as  such.  In 
view  of  their  strong  dependence  on  local  conditions 
they  are  not  included  in  the  present  modeling  effort. 

The  proposed  model  then  is  a  power  spectrum  of 
rotational  deflections  consisting  of  an  underlying 
random  walk  process  which  is  observable  over  ten 
decades  of  frequency,  superimposed  to  which  is  the 
spectrum  of  diurnal  random  fluctuations  as  described 
in  equation  (3). 

III.  Limitations  and  Uncertainties 
of  the  Model 

The  technique  of  power  spectral  representation 
of  time  series  has  certain  limitations  which  must  be 
kept  in  mind  in  interpreting  the  data.  First,  it  is 
valid  only  for  stationary  series,  the  statistics  of 
which  do  not  change  with  time.  Thus,  care  must  be 
taken  to  avoid  drawing  conclusions  from  data  which 
were  measured  during  periods  when  unusual  distur¬ 
bances  such  as  earthquakes  or  similar  events  take 
place.  However,  practical  experience  indicates  that 
periods  with  evidently  stationary  characteristics 
prevail  at  test  sites  so  that  a  meaningful  analysis 
can  be  conducted  with  reasonable  care.  Furthermore, 
it  must  be  recalled  that  power  spectral  densities 
computed  from  a  finite  amount  of  data  are  only 
estimates  of  the  true  values.  The  uncertainty  of 
these  estimates  can  be  as  high  as  ±100%  if  only  one 
set  of  measurements  is  available.  A  detailed 
analysis  of  these  uncertainties  is  given  in  Refer¬ 
ence  1.  In  general,  the  uncertainty  in  a  calcu¬ 
lated  spectrum  increases  with  decreasing  frequency 
which  may  lead  to  an  erroneous  determination  of  the 
slope  of  the  spectrum.  The  technique  applied  in 
this  paper,  namely,  to  line  up  spectra  measured  in 
different  frequency  intervals  with  different 
instrumentation  tends  to  de-emphasize  the  slope 
uncertainties  in  the  individual  spectra  and  to 
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provide  information  on  the  wide-range  slope  trends 
of  the  phenomenon  under  study. 

In  addition  to  the  statistical  uncertainties 
discussed  above,  the  inherent  uncertainties  of  the 
measurements  must  be  taken  into  consideration.  At 
the  low  end  of  the  spectrum,  the  uncertainties  con¬ 
nected  with  the  azimuth  calibration  of  a  reference 
monument  from  Polaris  observations  play  an  important 
role.  While  the  internal  random  noise  of  most 
tiltmeters  is  sufficiently  low  for  the  purposes 
discussed  here,  many  of  them  show  a  significant 
temperature  sensitivity  and  extreme  care  must  be 
taken  to  avoid  erroneous  results  due  to  this  error 
source. 

The  claim  of  the  proposed  model  to  be  valid 
within  the  stated  tolerance  limits  for  representa¬ 
tive  test  sites  seems  to  be  contradicted  by  a  mea- 
surem'ent  recently  made  at  the  Agassiz  test  station 
in  Harvard,  MA  by  the  C.  S.  Drapej  Laboratory.  The 
power  spectrum  of  these  measurements  which  was  made 
available  to  the  author  through  private  communica¬ 
tion  by  Mr.  Koenigberg,  shows  values  which  are  near¬ 
ly  three  orders  of  magnitude  smaller  than  expected 
from  the  model  in  the  frequency  range  from  3  x  10-3 
Hz  to  3  x  10"2  Hz.  This  is  the  only  significant 
deviation  from  the  proposed  model  which  is  presently 
known  to  the  author.  However,  the  uncertainties  of 
the  spectrum  are  particularly  large  in  this  fre¬ 
quency  range  as  will  be  noted  by  inspection  of 
Figure  5.  This  is  the  transition  range  where  motions 
generated  by  quite  different  mechanisms  overlap  and 
thus  such  large  uncertainties  may  be  explainable. 

IV.  Applications  of  the  Model 

The  availability  of  such  a  quantitative  ana¬ 
lytical  model  is  a  valuable  tool  for  the  planning 
and  evaluation  of  gyro  tests  in  many  respects.  One 
important  application  is  that  it  provides  an  under¬ 
standing  of  the  growth  of  the  uncertainty  in  the 
knowledge  of  the  orientation  of  an  azimuth  reference 
monument  as  a  function  of  * ime  after  calibration. 

This  information  can  be  derived  from  ..:i  integration 
of  the  power  spectrum  over  the  frequency  interval 
from  +  *  to  1/T  where  T  is  the  time  passed  since 
the  last  calibration.  This  translation  between  the 
frequency  and  the  time  domain  is  based  on  the  prop¬ 
erty  of  the  Fourier  transform  that  the  lowest  fre¬ 
quency  which  is  theoretically  defined  in  an  obser¬ 
vation  interval  T  is  1/T.  Curve  A  ir.  Figure  7  shows 
the  uncertainty  variance  of  the  random  wall  compon¬ 
ent  which  increases  proportional  to  time.  Curve  B 
is  the  variance  connected  with  the  diurnal  peak  of 
the  model  which  dominates  the  uncertainty  growth 
process  within  the  time  interval  of  piactical  inter¬ 
est.  The  sharp  increase  of  the  uncertainty  in  the 
time  interval  up  to  one  lay  must  h-  need  The 
understanding  of  the  statistical  nature  of  the 
monument  motions  provided  by  the  model  together  with 
the  knowledge  of  the  uncertainties  of  Polaris  obser¬ 
vations  can  be  used  to  optimize  the  calibration 
process  of  the  monument . 

Another  important  application  of  the  model  is 
its  suitability  to  predi^'  t  -  output  of  a  gyro 
under  test  which  is  due  to  test  pad  motions.  Curve 
A  in  Figure  8  shows  the  resulting  power  spectrum  of 
the  earth  rate  components  wh.ch  thr  gyro  "pic«r.  up" 
due  to  its  deflections  from  a  neutral  tesi  orienta¬ 
tion.  Curve  B  shows  the  gower  spectrum  of  the  pad 
motion  angular  rates  which  the  input  -xi  >'  h. 
gyro  senses.  Curve  r  is  ».  ;  the  drift 


rates  due  to  the  pad  motion  angular  accelerations 
sensed  by  the  output  axis  of  a  representative  gyro. 
This  externally  induced  gyro  output  noise  spectrum 
must  be  separated  from  the  internal  gyro  noise,  the 
quantity  of  interest  to  be  determined  in  the  test. 
Independent  on  how  this  separation  of  the  noise 
sources  is  approached,  either  servo  stabilization 
of  the  test  platform  or  compensation  in  the  data 
evaluation,  the  availability  of  a  quantitative 
estimate  of  the  external  noise  statistics  is  an 
important  tool  in  the  optimization  of  the  respective 
process . 

Furthermore,  the  knowledge  of  this  external 
noise  spectrum  can  be  used  to  assess  the  deteriora¬ 
tion  of  other  tests  due  to  the  motions  of  the  test 
pad.  In  Reference  5  this  is  done  for  a  standard- 
torque-to-balance  test.  Studies  concerning  the 
effects  on  a  four-position  gyro  compass  process  are 
under  way. 

V.  Conclusions 

The  paper  demonstrates  that  the  attempt  to 
formulate  a  model  of  angular  test  pad  motions  of 
general  validity  is  not  unreasonable.  However, 
more  observations  are  required  to  obtain  better 
estimates  of  the  tolerances  of  the  model  and  its 
structural  details.  In  particular,  more  information 
is  needed  to  explore  the  nature  of  the  microseism 
peak  of  the  spectrum.  To  accomplish  this  it  is  pro¬ 
posed  that  relevant  data  which  are  either  available 
now  at  various  test  organizations  or  which  will  be 
collected  in  the  future  be  made  available  to  the 
Geokinetics  Sub-Committee  of  the  AIAA  Technical 
Committee  on  Guidance  and  Control.  In  this  way  a 
basis  could  be  established  for  the  preparation  of  a 
Summary  Paper  or  a  Progress  Volume  and  the  formula¬ 
tion  of  a  generally  recognized  model. 
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FIG.  Z  SPECTRA  FROM  Ilf4  -  1  Hz 

0  EW  Tilt  Bedford  Flight  Fac.,  CSDL;  Apr  73 
Hughes  TM  3  Tiltmeter  (Reference  6) 
O  EW  Tiltf  Holloman  AFB,  CIL  Bldg. ,  August  71 
X  NS  Tiltj  Geotech  Model  18279  Tiltmeters 
X  EN  Tilt  Holloman  AFB,  CIL  Bldg.  April  74 
ADL  Mod  III  Tiltmeter 


•  EW  Tilt  Bedford  flight  Fac.,  CSDL; 
l  Feb  -  Apr  73 

O  NS  Tilt  I  Hughes  TM  3  Tiltraeter 
J  (Reference  6) 

X  Composite  EW  -  NS  Tilt  Vector 

Holloman  AFB,  Lab  Bldg.  , 

June  67 

Geotech  Model  18279  Tiltmeter 
(Reference  7) 


I  ■  ■  I  Average  Azimuth  Fluctuations 

Anaheim,  California 

Autonetics  R  6  D  Azimuth  Reference 

Calibrations.  (Reference  9) 

A  May  -  Nov  67 

B  Dec  67  -  Feb  70 


FIG.  5  COMPOSITE  SPECTRUM 
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FIG.  6  PROPERTIES  Of  THE  MICROSEISMIC  PEAK 

:  Squared  Coherence  Function  of  Measure¬ 
ments  from  Two  Seismometers  Placed  at  a 
Distance  of  153  m. 

(Transcribed  from  Reference  3,  Figure  4) 
:  Ratio  of  Power  Spectral  Densities  of 
Translational  Acceleration  Measured  at  a 
Depth  of  668  m  and  on  the  Surface. 
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FIG.  7 

GROWTH  OF  PAD  ORIENTATION  UNCERTAINTY  WITH  TIME 


A  :  Influence  of  Random  Walk  Process 
B  .  Influence  of  Diurnal  Fluctuations 
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